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Abstract 

We study the current and the associated noise for the transport through a two-site 
molecule driven by an external oscillating field. Within a high-frequency approxi- 
mation, the time-dependent Hamiltonian is mapped to a static one with effective 
parameters that depend on the driving amplitude and frequency. This analysis al- 
lows an intuitive physical picture explaining the nontrivial structure found in the 
noise properties as a function of the driving amplitude. The presence of dips in the 
Fano factor permits a control of the noise level by means of an appropriate external 
driving. 
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1 Introduction 



The development of techniques to contact molecules and to drive a current 
through them [1,2,3] has opened the new and promising field of molecular elec- 
tronics [4,5]. In order to construct useful devices, however, it is not sufficient 
to have a current flowing through a molecule but one also needs to have the 
ability to control this current. This can be achieved in principle by means of 
the so-called single electron transistor setup where a gate electrode is placed 
close to the molecule. Applying a gate voltage then allows to influence the 
current across the molecule. In more complex circuits, the need for a large 
number of contacts or electrodes close to the molecule may constitute a major 
obstacle. In fact, already the implementation of a single gate electrode which 
creates a sufficiently strong field at the molecule is a demanding task [6,7,8]. 
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Therefore, other means of controlling the current through a nanosystem should 
be explored. One possibility is to replace the static field of a gate electrode 
by a suitable external ac field. Recent theoretical work [9] has demonstrated 
that, by using a coherent monochromatic field, one should indeed be able to 
control the electrical current flowing through a nanosystem connected to sev- 
eral leads. In an extension [10] of this work, it was demonstrated that even 
the noise level can be suppressed by an appropriate driving field. 



In the present paper, we represent the molecule by a two-site system under 
the influence of an external high-frequency field and coupled to leads. Such 
a model is not limited to describe electrical transport through molecules but 
may also be applied to other situations like coherently coupled quantum dot 
systems [11] irradiated by microwaves. 



In Ref. [10], a Floquet approach was employed to derive exact expressions for 
both, the current and the associated noise for the transport through a non- 
interacting nanosystem in the presence of an arbitrary time-periodic field. A 
study of the Fano factor, i.e. the ratio between noise and current, revealed a 
suppression at certain values of driving amplitude and frequency. To achieve a 
better physical understanding of this phenomenon, we here consider the prob- 
lem within the high-frequency regime, which allows us to approximate the 
driven system by a static one with renormalised parameters. The structure 
observed for the Fano factor can then be understood in terms of three differ- 
ent scenarios. For small effective intramolecular hopping matrix elements, the 
system itself acts as a bottleneck, while in the opposite limit, the two con- 
tacts form a two-barrier setup. In between, when the hopping matrix element 
is of the order of the system-lead coupling strength, the barriers effectively 
disappear, leading to a suppression of the Fano factor. 



In the next section, we introduce our model, consisting of two sites subject 
to an external oscillating field and coupled to two leads. We then discuss this 
model in the static case, deriving explicit expressions for both the current and 
its noise. In Sect. 4, expressions for the effective hopping matrix element and 
the electron distribution functions in the leads are determined within the high- 
frequency approximation. These results are used in Sect. 5 to compute current, 
noise and the corresponding Fano factor. A comparison with results based 
on the exact expressions of Refs. [10] demonstrates the validity of the high- 
frequency approximation for not too low frequencies. This allows a physical 
interpretation of the observed features in terms of a static model. 
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Fig. 1. Level structure of the nanoscale conductor with two sites. Each site is coupled 
to the respective lead with chemical potentials hl and \ir = hl + eV. 



2 The model 



In the following, we consider the setup depicted in Fig. 1, which we describe 
by the time-dependent Hamiltonian 

H{t) = H sys t em {t) + H \ ea d s + -^contacts • (1) 

The first term on the right-hand side, 

H systcm (t) = -A(c\c 2 + c\ci) + ^(c\ci - c\c 2 ) cos(fit), (2) 

represents the driven two-site system, where electron-electron and electron- 
phonon interactions have been disregarded. The fermion operators c n and 
c4, n — 1,2, annihilate and create, respectively, an electron at site n. Both 
sites are coupled by a hopping matrix element A. The applied ac field with 
frequency Q = 2ir/T results in a dipole force given by the second term in the 
Hamiltonian (2). The amplitude A is proportional to the component of the 
electric field strength parallel to the system axis. 

The electrons in the leads are described by the Hamiltonian 

-Pleads = XX e ^<7 C Lq C Lq + € Rq C Rq C Rq)i ( 3 ) 
9 

where (c^) creates a spinless electron in the left (right) lead with momen- 
tum q. The electron distribution in the leads is assumed to be grand canonical 
with inverse temperature (3 = 1/ksT and electro-chemical potential Hl/r- An 
applied voltage V corresponds to Hr — Hl = eV, where — e is the electron 
charge. 
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The tunnelling Hamiltonian 

-^contacts = E (^ig C L<j C l + V Rq c ] Rq c 2 ) + h.c. (4) 
i 

establishes the contact between the sites and the leads, as sketched in Fig. 1. 
The system-lead coupling is specified by the spectral density 

F e (E) =2tt^|^| 2 5(E -e qt ) (5) 
q 

with £ — L,R. Below, we shall assume within a so-called wide-band limit that 
these spectral densities are energy independent, T^E) = IV 



3 Transport through a static two-site system 



We start by deriving expressions for current and noise for a static two-site 
system coupled to two leads, setting A = in the Hamiltonian (2). Solving 
the Heisenberg equations of motion for the lead operators, we obtain 



c Lq (t) = cz^oK^^'^ -^Jdt' 

to 



-ie Lq (t-t')/h 



Cl (0 



(6) 



and a corresponding expression for CR q {t) with L replaced by R and C\ by c 2 . 
Inserting (6) into the Heisenberg equations of motion of the two-site system 
and exploiting the wide-band limit, one arrives at 



r 

C! 4 Ac2 -2t Cl+ ^ (t) ' 

c 2 4 Aci -^f C2+ ^ (t) - 



(7) 
(8) 



For a grand canonical ensemble, the operator-valued Gaussian noise 

= ~\ E K ex P \-T e ^ - *o)l C M (9) 



obeys 



Ut)) = o, 



(10) 

fill 



where /^(e) = {1 + exp[/5(e — /x^)]} 1 denotes the Fermi function with chemi- 
cal potential £ — L,R. In the asymptotic limit to - ► — oo, the solutions of 
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Eqs. (7) and (8) read with n = 1, 2: 

oo 

c n (t) = Jdr {G u1 (t) £ L (t - r) + G n2 (r) £ R (t - r)}. 



(12) 



In the wide-band limit and for equal system-lead coupling, T e = T, the prop- 
agator is given by 



G(r) = e- r ^ 2 



^cos(At) isin(Ar) 
Jsin(Ar) cos(At) 

where Q(r) is the Heaviside step function. 



Sir) 



(13) 



The operators corresponding to the currents across the contacts £ = L, R are 
given by the negative time derivative of the electron number Ng = J2 q c l q c e q 
in the respective lead, multiplied by the electron charge — e. For the current 
through the left contact one finds 



^r L c\(t) Cl (t) - e{c\(t)t L (t) + d(*)ci(*)} 



(14) 



with a corresponding expression for Iii(t). In the stationary limit, to ~~ *■ — oo, 
the mean values of the currents across the two contacts agree and we obtain 



(h) 



2nh 



dE 



f R (E)-f L (E) T(E). 



(15) 



In the wide-band limit, the transmission T(E) can be expressed in terms of 
G 12(E), i.e. the Fourier transform of the propagator G 12 (t), as 

T(E) = T L T R \G 12 (E)\ 2 . (16) 

Making use of the propagator (13), the transmission for T e = T becomes 

T 2 A 2 



T(E) 



|(£-ir/2) 2 - A 2 | 2 ' 



(17) 



The noise of the current through contact i is given by the symmetric auto- 
correlation function of the current fluctuation operator AIg(t) = Ig(t) — {Ig(t)). 
It is possible to characterise the noise strength by its zero frequency component 

S = \ [ + °° dt(AI e (t)AIg(0) + A/,(0)A/,(t)>, (18) 
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which is independent of the contact £. The quantity S may be expressed in 
terms of the transmission function T(E) as [12] 



2nh 



J dE {T(E)[f L (E)[l - f L (E)} + f R (E)[l - f R (E)\ 
+ T(E) [l - T(E)} [f R (E) - f L (E)fy 



(19) 



Two contributions to the zero-frequency noise S have to be distinguished: 
The first term is a temperature-dependent equilibrium noise according to the 
dissipation-fluctuation theorem [13] and dominates for (3eV <C 1. In contrast, 
for large voltages fieV ^> 1, the main contribution to the noise stems from the 
second term. This so-called shot noise has its physical origin in the discreteness 
of the charge carriers. 

We now consider voltages larger than all other energy scales in the problem. 
As a consequence, the current noise will entirely be due to shot noise. Fur- 
thermore, in this limit, the results for current and noise will not depend on 
temperature. In the expression (15) for the current, the difference of the Fermi 
distributions then practically equals one for energies where the transmission 
is nonvanishing. The current thus reads 



A 2 



27rft 2h A 2 + (r/2) 



2 ' 



(20) 



where T = J dET(E) is the total transmission. With the same argument we 
find from (19) for the current noise 



5*00 — 



e 2 r 2A 2 (r 4 - 2r 2 A 2 + 8A 4 ) 



h 



(4A 2 + r 2 y 



(21) 



The relative noise strength can be characterised by the so-called Fano factor 
F = S/el which, in the infinite voltage limit, becomes 



r 4 - 2r 2 A 2 + 8A 
(4A 2 + r 2 ) 2 



(22) 



In Fig. 2, the Fano factor F^ is depicted as a function of the ratio of the 
tunnelling matrix element A and the level width T. For weak system-lead 
coupling r < A, the two contacts between the two-site system and the leads 
form the limiting step of the transport process. We effectively arrive at trans- 
port through a double-barrier system with a Fano factor Foq = 1/2 [14]. On 
the other hand, for T ^> A the two sites hybridise with the adjacent lead and 
effectively only a single barrier remains. This yields a Fano factor F^ = 1. At 
the crossover between these two regimes, the channel is optimally transparent 
and, consequently, the Fano factor assumes a minimum. From the expression 
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I 1 1 1 1 1 

\pk 1 2 3 4 

A/r 

Fig. 2. Fano factor = Soo/el^ as a function of A/r. For A <C T, the bottleneck 
of the transport is the tunnelling process between the two sites yielding a Fano factor 
Foo = 1. In the opposite limit A > T, we obtain transport through a double-barrier 
structure with a corresponding Fano factor F^ = 1/2. In the intermediate regime 
the Fano factor assumes a minimum at the position indicated in the plot. 

(22), we find the optimal hopping matrix element A = ^5/12 T yielding a 
minimal Fano factor of = 7/32. We remark that the minimum decreases 
further if the number of sites in the system is increased [10]. 



4 High-frequency approximation 



Let us now turn back to the original time-dependent problem. We will compute 
within a high-frequency approximation [15] the current through this driven 
system and the corresponding current noise. Results valid for arbitrary driving 
amplitudes A can be obtained by the following procedure which is justified in 
the Appendix on the basis of Floquet theory. 

First, we introduce the interaction picture with respect to the driving which 
for the problem at hand is obtained by means of the unitary transformation 

U (t) = exp (-i— (cjci - c\c 2 ) sin(m)) • (23) 

This yields the new system operators 

ci, 2 (0 = U^(t)c h2 U (t) = c 1>2 exp (V ^ sin(QO) , (24) 

where the upper sign corresponds to site 1. To a good approximation, the 
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dynamics can then be described by the time-averaged system Hamiltonian 

^system = -A cff (c t 1 C 2 + c|ci). (25) 

Thus, within a high-frequency approximation, the driven two-site system acts 
as a static system with the effective hopping matrix element 

A cfr = J (A/hQ)A, (26) 

where Jo is the zeroth order Bessel function of the first kind. The driving 
amplitude A and frequency Q can now be chosen such that A e g vanishes 
and consequently tunnelling between the two central sites then no longer oc- 
curs [16,17,15]. 

Proceeding as in Sect. 3, the influence of the leads after the transformation 
(23) can be described by fluctuation operators. For the left lead one finds 

VL(t) = ~ E ^L>xp [-i (e Lq {t - t ) + A sin(Qt))] c Lq (t ) (27) 
with the correlation function 



'"' k,k' 

xJ k (A/2hQ)J k ,{A/2hQ) e"^*-*')"* (28) 

and corresponding expressions for the right lead. Since we are interested in 
the average current and the zero-frequency noise, i.e. low-frequency transport 
properties, we can neglect the T-periodic contribution to the correlation func- 
tion (28) and, thus, average its t-dependence over the driving period. Then, 
the correlation function assumes the form (11) like in the static case but with 
the Fermi function replaced by the effective distribution function 

fi, cS (E)= £ J 2 k (A/2hQ)ME + khQ). (29) 

k=— oo 

The different terms in this sum describe processes where an electron of en- 
ergy E is transferred from lead I to the system under absorption (emission) 
of | A; | photons for k < (k > 0). These processes are weighted by the square 
of the kih order Bessel function of the first kind. 

Having approximated the originally time- dependent problem by a static one 
with an effective hopping matrix element and an effective distribution function, 
we can calculate the transmission, the current, and the zero frequency noise of 
the driven system with the formulae derived in Sect. 3 for a static situation. 
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Fig. 3. Typical energy dependence of transmission T{E) (solid line) and effective 
distribution function fe, e fi{E) (dashed line) which allows to replace the distribution 
function by (30) in the expressions for current and noise. 

In the limit of very large voltages and for energies where the transmission 
is nonvanishing, the effective distribution functions in the left and right lead 
become again zero or one, respectively. As a consequence, the time-averaged 
current and the zero-frequency noise are given by the expressions (20) and 
(21) with the replacement A — > A e g. We denote the current and the noise 
in this limit by 1^ and Soo, respectively. As pointed out above, there exist 
driving parameters where the effective hopping matrix element (26) vanishes. 
As a consequence, no current can flow through the system under these cir- 
cumstances [9,10]. 

The case of a finite voltage requires a more detailed inspection of the distri- 
bution functions fe fiS and the effective transmission T(E) sketched in Fig. 3. 
In the high-frequency regime under study here, the width W of the transmis- 
sion function is much smaller than Ml. Furthermore, the effective distribution 
function /^ e fr is nearly constant for energies E separated by at least 1/(5 from 
the steps at E — + kMl. Therefore, unless a step in /^ e ff occurs close to 
E = 0, the effective distribution functions in the current and noise expressions 
can be replaced by their value at E — 0, i.e. 



(30) 



Thus, the time-averaged current and the zero-frequency noise are given by 

(31) 
(32) 

where K(V) denotes the largest integer not exceeding eV/2Ml. Note that the 



/ K(y) \ 

I = lJj 2 Q (A/2hQ) + 2 ]T J 2 k (A/2hQ) , 

^ k=l ' 

s = + [Jl(A/2m) + 2 p J 2 k (A/2hn)j (X - 
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Fano factor F = S /el for fixed A/Q reaches its minimal value in the infinite 
voltage limit. Since Jk{x) ~ for x > k and J2k Jl( x ) = 1> the dc current and 
the zero frequency noise are well approximated for A < eV by their asymptotic 
values for infinite voltage, I as 1^ and S rs S^. We remark that, in contrast 
to the static case, the result (32) contains contributions stemming from the 
first term in the noise expression (19) even in the zero-temperature limit. 



5 Comparison with exact results 

Figures 4b-d depict by solid lines the time- averaged current, the zero-frequency 
noise and the Fano factor at zero temperature obtained numerically within the 
Floquet approach of Ref. [10] for the relatively large voltage V = 48A/e. This 
particular value of the voltage has been selected to avoid the chemical poten- 
tials to lie close to multiples of Kl. A comparison of these numerically exact 
results for current and noise with the approximate expressions (31) and (32) 
depicted by dashed lines shows a good agreement for the parameters chosen. 
The agreement improves with increasing frequency: already for Q = lOA/ft, 
it is found that the exact and approximate results can practically no longer 
be distinguished. 

The exact numerical results show strong suppressions of both, the current and 
the noise for certain driving amplitudes. This behavior can be explained within 
the high-frequency approximation presented in Section 4: Whenever the ratio 
A/Ml corresponds to a zero of the Bessel function J , the effective hopping 
matrix element A e g vanishes (cf. Fig. 4a) and consequently the current and 
the noise become zero. Note that the exact result exhibits still a residual 
current and noise. The suppressions of the current and noise lead to peaks of 
the Fano factor F. For sufficiently small driving amplitudes, these peaks are 
accompanied by minima which correspond to |A e ff| ~ \J~hjv2T indicated by 
the dotted line in Fig. 4a. 

For driving amplitudes A < eV, the finite voltage results (31) and (32) for 
the current / and the noise S are well described by the results (20) and (21) 
for infinite voltage with A replaced by A eff . In this regime, the Fano factor 
reaches maxima F = 1. In contrast, for larger driving amplitudes A > eV, we 
find a Fano factor larger than that predicted by (22), as discussed below (32). 
In particular, the Fano factor can assume values F > 1. 

Finally, we consider in Fig. 5 the case of intermediate voltages such that 
A, T < eV < 2hQ. Then, only the zero photon channel contributes and hence 
the current I = I 00 Jq{A/2MT) is considerably lower than for large voltages. 
Now, in addition, a new type of suppression appears at twice the amplitude 
compared to the suppressions discussed above. The physical reason for this new 
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Fig. 4. (a) Effective hopping matrix element |A e g|, (b) time-averaged current I, 
(c) zero- frequency noise S, and (d) Fano factor F = S/el as a function of the 
driving amplitude A. Shown are the numerically exact results (solid lines), the 
approximative results (31) and (32) for finite voltage (dashed lines), and the infinite 
voltage results (20) and (21) with A replaced by A c g- (dotted lines). The coupling 
strength is T = 0.5 A, the driving frequency is Q = 5A/h, and the voltage reads 
V = 48A/e. The dotted line in (a) marks the value V for which the Fano 

factor assumes its minimum. 
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Fig. 5. Time-averaged current I as a function of the driving amplitude A. Shown 
are the numerically exact result (solid line), the approximate result (31) for finite 
voltage (dashed line) and the infinite voltage result (dotted line). The coupling 
strength is T = 0.5 A, the driving frequency is Q = 5 A/ft, and the voltage reads 
V = 5A/e. 

kind of suppression lies in the fact that the effective distribution functions in 
the two leads are equal at the relevant energies and therefore no dc current 
can flow. Nevertheless, the noise remains finite and, consequently, the Fano 
factor diverges. 



6 Conclusions 

We have presented a high-frequency approximation for the charge transport 
through a driven two-site system. Within this scheme, the time-dependent 
Hamiltonian and the lead correlation functions are transformed to an appro- 
priate interaction picture and subsequently time-averaged over the driving 
period. In the resulting equations, the hopping matrix element of the two-site 
system and the electron distributions of the attached leads are replaced by 
effective ones which depend on the driving parameters. 

This static picture allows to gain profound physical insight into the structure 
present in current and noise as a function of the driving parameters. For small 
effective hopping matrix elements, the barrier between the two sites of the 
system dominates, leading to shot noise with F = 1. In the opposite limit, the 
contacts form a double-barrier system corresponding to a Fano factor F — 1/2. 
Between these two situations, effectively no barrier exists in the transport path 
and the Fano factor is further reduced to a minimal value of F — 7/32. 

The results of this work demonstrate that the control of a current through a 
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molecule by means of a time-periodic driving provides a viable alternative to 
the traditional single electron transistor setup based on a gate electrode. Our 
approach allows to minimise the number of electrodes close to the molecule. 
Furthermore, a suitable choice of the transistor's working point permits to 
operate in a low noise regime with a small Fano factor. These features inherent 
to the driven setup may prove useful for the development of novel molecular 
electronics devices. 
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A Driven quantum systems in high-frequency approximation 

In this Appendix, we review a common perturbative approach for the treat- 
ment of periodically time-dependent quantum systems and thereby justify the 
high-frequency approximation employed in Section 4. 

A standard technique for the study of periodically time-dependent Hamil- 
tonians H(t) = H{t + T) is the so-called Floquet approach [18,19,20]. It 
starts out from the fact that a complete set of solutions of the corresponding 
Schrodinger equation is of the form \ip a (t)) = e~ ltat ^ h \(j) a (t)) where the Flo- 
quet states \4> a (t)) = \<p a (t + T)) obey the time-periodicity of the Hamiltonian. 
The Floquet states and the quasi-energies e a are eigenstates and eigenvalues, 
respectively, of the Hermitian operator TL = H(t) — ihd/dt defined in a Hilbert 
space extended by a periodic time coordinate. We emphasise that already the 
Floquet states from a single Brillouin zone —Ml/ 2 < e a < hfl/2 form a com- 
plete set of solutions. A Floquet ansatz essentially maps the time-dependent 
problem to an eigenvalue problem and, thus, it is possible to employ all ap- 
proximation schemes known from time-independent quantum mechanics, in 
particular perturbative schemes for the computation of eigenstates. 

Here, we consider the special case of a time-dependent Hamiltonian of the 
form 



where f(t) is a T-periodic function with zero time-average. If HQ is much larger 
than all energy differences in the spectrum of H 1 , the following Schrodinger 



H = H f(t) + H 



i 
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perturbation theory can be employed for the computation of the Floquet states 
[18,19]: It is assumed that for H Q all eigenstates \tp a ) and eigenenergies E a are 
known. Then, the unperturbed Floquet Hamiltonian H, = H f(t) — ihd/dt 
has the eigensolutions 

\<j) k a (t)} = exp (- l -E a F{t) + ifcftf) \<p a ) (A.2) 

with eigenvalue khQ. Here, dF(t)/dt = fit) and k is an arbitrary integer. 
Note that F(t) satisfies the T-periodicity of the field since the time average of 
fit) vanishes. Thus, k defines a degenerate subspace of the extended Hilbert 
space. In each degenerate subspace, the matrix elements of the perturbation 
read 

{H,U^^[ dt(<t>l{t)\HM k M. (A.3) 

Therefore, to first order in Hi /Ml, the Floquet states and the quasienergies for 
the Hamiltonian (A.l) are obtained by diagonalising the perturbation matrix 
(A.3). 

Following (A.2), the basis states \ip a ) and |<^(t)) are related by the unitary 
transformation 

U (t) = exp (-±.H F(tj) (A.4) 

as exp (— ikflt) \(f) k a {t)) = U (t)\ip a ). For the Hamiltonian (2), (A.4) corresponds 
to the unitary transformation (23). 

Within the regime of validity of the perturbative approach, the problem is 
therefore described by the static Hamiltonian 

#i = ^ C dtUl{t)E x U Q {t). (A.5) 
1 Jo 

Note that after the transformation with (A.4), the amplitude of the oscillating 
part of the new Hamiltonian ul(t)H-JJo(t) is no longer governed by H , but 
rather by H 1 . Thus, a perturbative treatment of the oscillating part is (almost) 
independent of the original driving amplitude in H . 

A particular example for a high-frequency approach along these lines is a 
particle moving in a one dimensional continuous potential under the influence 
of a dipole field, i.e. Hi = p 2 /2m + V(x) and H = fix. Then, (A.4) constitutes 
a gauge transformation and results in a Hamiltonian which is again of the 
form (A.l). A second transformation of the type (A.4) yields a periodically 
accelerated potential and defines the so-called Kramers-Henneberger frame 
[21,22]. 
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